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ABSTRACT 

Context. The complexity of the ISM is such that it is unlikely that star formation is initiated in the same way in all 
molecular clouds. While some clouds seem to collapse on their own, others may be triggered by an external event 
such as a cloud/flow collision forming a gravitationally unstable enhanced density layer. 

Aims. This work tests cloud-cloud collisions as the triggering mechanism for star formation in the Serpens Main 
Cluster as has been suggested by previous work. 

Methods. A set of smoothed particle hydrodynamics (SPH) simulations of the collision between two cylindrical 
clouds are performed and compared to (sub)millimetre observations of the Serpens Main Cluster. 
Results. A configuration has been found which reproduces many of the observed characteristics of Serpens, 
including some of the main features of the peculiar velocity field. The evolution of the velocity with position 
throughout the model is similar to that observed and the column density and masses within the modeled cloud 
agree with those measured for the SE sub-cluster. Furthermore, our results also show that an asymmetric collision 
provides the ingredients to reproduce lower density filaments perpendicular to the main structure, similar to those 
observed. In this scenario, the formation of the NW sub-cluster of Serpens can be reproduced only if there is a 
pre-existing marginally gravitationally unstable region at the time the collision occurs. 

Conclusions. This work supports the interpretation that a collision between two clouds may have been the trigger 
of the most recent burst of star formation in Serpens. It not only explains the complicated velocity structure seen in 
the region, but also the temperature diff'erences between the north (in "isolated" collapse) and the south (resulting 
from the shock between the clouds). In addition it provides an explanation for the sources in the south having a 
larger spread in age than those in the north. 

Key words. Stars: formation - Stars: individual: Serpens - ISM: clouds - ISM: kinematics and dynamics - 
Methods: numerical 



1. Introduction 



Observations of molecular clouds indicate that stars 



form in dense, clump y filaments (e.g. Andre et al. 2010 
Molinari et al.||2010 ). However, the processes control- 



ling the formation of these structures and their role in 
triggering specific star formation episodes are widely 
debated. One process could simply be gravity con- 
trolled: quasi-static or dynamic gravitational collapse 
of non-spherical clouds can easily generate filamen- 
tary struc tures which subsequently fragment to form 
stars (e.g. |Palla & Stahler„2002, ,Hartmann & Burkertj 
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|2007| ). The contraction of an idealised isolated turbulent 
cloud leads to the formation of small scale filaments 
(~ 0. 1 pc long) within which stars begin to form, such as 



Klessen 



shown by simulations from |Bate et aL](|2003| ) 
[eTaLl (|2n05]); |Batel (|2009a| and|2009b)). Alternatively, 
external triggers can compress molecular gas to gen- 
erate high-density fragmented filaments. Such triggers 
can be of very diflTerent natures: ionisation/shock fronts 
around OB stars (e.g. Elmegree n & Lada]|1977[ |Dale| 



et al. [2007 ), large scale colliding flows generated by 



supemovae/galactic shear (e.g. [Heitsch et al. 2008 ), or 
even molecular cloud collisions (e.g. [Gittins et al.,2003] 
Kitsionas & Whitworth] [20071 |Anathpindika||2009a[ 



2009b| ). 



Each of these processes have been tested using nu- 
merical simulations, and their respective eflftciency in 
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forming stars have been investigated. For some local 
clouds, we are beginning to acquire sufficiently de- 
tailed observations to generate a more informed picture 
of how star formation has progressed in those clouds. 
With these observations it becomes possible to per- 
form dedicated simulations to model and better un- 
derstand the processes at work. Only few such studies 
have been per formed so far. For instance, [ Ballesteros^ 
Paredes et al.|(l99 9 ) used the results of their turbulent 
colliding flow simulations to argue that Tauru s-Auriga 
formed via such a process; Peretto et al. ( 2Q07| ) used hy- 
drodynamical simulations to model the global collapse 
and fragmentation of the NGC 2264-C protocluster; 
and Heits ch et al.| ( |2Q09| compared flow-driven models 
for the formation of isolated molecular clouds with the 
Pipe nebula. In this paper, we use numerical simulations 
to investigate whether the recent star formation episode 
observed in the Serpens Main Cluster could have been 
triggered b y the collision of two mole cular clouds, as 
proposed in Duarte-Cabral et al. ( |201Q ). 

Serpens is an interesting region in several aspects. 
It is a filamentary region with two compact young star 
forming clumps. The gas emission reveals an inter- 
esting velocity structure and has an unusual tempera- 
ture structure ( [Duarte-Cabral et al.|2010| ). The observa- 
tions point to remnant signatures of the trigger behind 
the current young on-going star formation of Serpens. 
Duarte-Cabral et aL] pOlOj suggested that this region 
appears to have been affected by a collision of two 
clouds or flows. Such a scenario is also suggested by 
recent studies of the global dynamic s of Serpens and 
other star forming regions (e.g. Hi guchi et al. 120 10} 
[Schneider et al.|2QlQt|GaIvan-Madrid et al.|20 1ol 

Here we use smoothed particle hydrodynamics 
(SPH) simulations of cloud-cloud collisions to attempt 
to reproduce the column density distribution and dy- 
namical properties of the Serpens Main Cluster. A brief 
summary of the observationally inferred characteris- 
tics of Serpens and the motivation for the simulations 
are presented in Section |2] Section [3] explains the SPH 
code, the initial conditions and the geometrical configu- 
rations adopted. The results of the calculations are pre- 
sented in Section]?] Section |5] presents a discussion of 
these results, where we compare the numerical results 
with the observations and choose the best- fit scenario to 
represent the region. The final remarks and conclusion 
are laid out in Section [6] Finally, in Appendix |A| (avail- 
able online only) we present temporal snapshots for all 
the models presented in this work. 



2. The structure of the Serpens Main Cluster 
region 

Serpens Main Cluster is a low mass star forming re - 

[19961 ). 



gion at 260 pc from the Sun (Straizys et al 



It comprises two compact protoclusters, hereafter re- 
ferred to as sub-clusters, lying in a 0.6 pc long fila- 
mentary structure, along a NW - SE direction. At first 
sight, the NW and SE sub-clusters appear rather sim- 
ilar. The dust emission ( [Davis et al.|[1999[ and Fig. 



NW 

sub-cluster 




600 

MJy/str 



SE 

^sub-cluster 



1 8^30'"00^ 50^ 
RA (J2000) 

Fig. 1. Serpens region as seen in 850 ji m dust continuum 
emission with SCUBA ( [Davis et al.|1999| ) in contours, at 0.4, 

0. 6. 1. 1.4. 1.8. 2.4 a nd 5 Jy beam~^ . In"g rey scale is the 24 yum 
from Spitzer MIPS ( [Harvey et al.|2007| ). All the sources seen 
in this image, both 24 fim and 850 fim sources, are Class and 

1, with only a few flat- spectrum sources. The 850 jim sources 
are labeled. 



[T]), shows that they have similar masses within simi- 
lar sized regions: ~ 30 M© in 0.025 pc^ eac h, as es- 



timated from both gas emission from C^^O ( Duarte 
[Cabral et al.[|2010] ) and 850 yum continuum emission 



from SCUBA. Furthermore, they both have sources at 
roughly the same stage of evolution. Class and Class I 



protostars, powering a number of outflows ( (Davis et al. 
[T9991 [Graves etal.|2010D . 



Previous studies (jKaas et aL]|2004[ [Harvey et al. 



'2006') had estimated an average age of 10^ yr for these 
sources, possibly marking a triggering event. Older 
sources are also found in the field, but they appear more 
dispersed, with no obvious connection to the current 
protostars, and no longer embedded nor surrounded by 
any cold dust seen at submillimetre wavelengths. It is 
thought these belong to a diflTerent burst of star forma- 
tion which occurred 2 Myr before the current one, and 
that these sources are now in the field as an unbound 
cluster. 

In contrast to the similarities of the dust emission, 
the velocity structure and molecular emission from each 
sub-cluster are strikingly diflTerent. The NW is a well 
"behaved" sub-cluster with one main velocity through- 
out, matching the systemic velocity of the cloud with 
only a smooth velocity shift towards the edges of the 
sub-cluster. This is shown in the C^^O J=1^0 position- 
velocity (PV) diagram of Fig. [2] (top) and on Fig. [3] 
which shows the structure of the higher and lower ve- 
locity components of Serpens. In this NW sub-cluster, 
the gas and dust emission have similar distributions, 
with the gas emission peaks corresponding to dust 
emission peaks. In addition, Spitzer 24 yum sources are 
coincident with many of the 850 yum sources, hereafter 
submillimetre sources (Fig.[T]). 
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Fig. 2. Position- velocity diagrams at constant declination 
of C^^O J=1^0 (grey scale and contours) as examples of 
the typical velocity structure in the NW sub-cluster (top) 
and the SE sub-cluster (bottom) ( Duarte-Cabral et al.|2010i ). 
Right Ascension varies from 18^30'"06^ to 18^29'^46^ for all 
panels. Declinations are 1° 16^48'' (top panel) and 1° 12^38'' 
(lower panel). The positions of these cuts in the map are 
shown as white lines in Fig. |3] The dashed lanes labeled as 
"dust" represent the regions whose S50/um emission is above 
0.6 Jy beam"^ 



In contrast, the SE sub-cluster is quite different and 
its velocity structure is more complex. Optically thin 
tracers show locations with a broad single component 
(at the north of the SE sub-cluster) which then splits 
into two clearly separated components (Fig.|3]and lower 
panel of Fig.|2]). At the southern end of the sub-cluster, 
the filamentary structure is then only seen in one ve- 
locity component al ong the line of sight (lower PV di- 
agram of Fig. 17 of Graves et al.||201Q| ). Furthermore, 
the gas and submillimetre emission peaks do not coin- 
cide. Rather the gas emissi on appears to peak be tween 
the submillimetre sources (Duarte- Cabral et aL][20 1Q). 
Finally, there are only a few 24 jjm sources associ- 
ated with the submillimetre dust continuum emission 
(Fig. [T]). The overall picture resembles that of a region 
where star formation is an on-going process, with some 
younger sources (the purely submillimetre sources) and 
others older (the purely 24 yum sources), unlike the NW 
sub-cluster where the sources appear to be all at the 
same evolutionary stage. 

A final diff'erence between the two Serpens sub- 
clusters is the gas temperature ( jDuarte-Cabral et al. 
|2010| ). The NW has a very uniform temperature around 
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Fig. 3. C^^O J=l-0 emission separated into blu e and red com- 
ponen ts, using the line fitting and splitting by Duarte-Cabral| 
|et al.| |2010 ). The blue emission corresponds to the total in- 
tegrated intensity of the low velocity cloud and the red emis- 
sion corresponds to the high velocity cloud. The submillime- 
tre sources from the 850 jim emission ( [Davis et al.|1 999 ) are 
marked with white triangles and serve as a guide to the loca- 
tion of the dust emission. The position of the PV diagrams in 
Fig. [2] are shown by the white solid lines, and the main axis 
of the blue shifted and red shifted cloud is shown as blue and 
red dashed line respectively. 



lOK, while the SE has both higher and more varied tem- 
peratures, ranging between 10 and 20K. In addition, the 
temperature peaks are not coincident with any submil- 
limetre source but rather between them, in the region 
where the two velocity components overlap. 

Thus the diff'erences between the two sub-clusters 
indicate that they have had a complicated history. A sce- 
nario capable of explaining the trigger of Serpens star 
formation has to reproduce the inhomogeneities in the 
sources' age, velocity structure and temperature distri- 
bution. The collision of two filament-like clouds, col- 
liding only over a portion of their length could provide 
such a trigger. 

3. SPH Simulations 

To test the proposal that the structure and star formation 
in the Serpens Main Cluster results from triggering by 
a cloud-cloud collision we have performed a set of SPH 
simulations which are compared in detail with millime- 
tre and submillimetre observations of Serpens. 

3.1. Numerical Code 

The calculations in this paper are performed using a 
SPH code based on a version by |Benz et al.| ( [199Q| ), but 
has been subje ct to substan tial modifications, including 
sink particles (|Bate[1995 ), variable smoothing lengths 
([Price & Monaghan[2005 J and magnetic fields ( Pric e &| 
Bate|2007| ). The code has been frequently used for sim 



ulations of star formation (e.g. [Dobbs et al.|2QQ6[ Bate 
[2009a| ). 
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We perform calculations which include the hydro- 
dynamics and self gravity, but do not include magnetic 
fields. We use sink particles, which are inserted in re- 
gions of high density (10"^^ g cm"^) that are undergo- 
ing collapse, to represent protostars. However for all the 
analysis in this paper we use the frame of the simulation 
when the first sink particle appears, i.e. before stellar 
feedback is likely to have an eff'ect. Therefore our re- 
sults are not dependent on the details, or dynamics of 
the sink particles. 

For all the calculations we adopt an isothermal 
equation of state. Again, we are only interested in the 
evolution of the clouds until star formation takes place. 
So this simplification is reasonable given that we con- 
sider the structure of the gas prior to stellar feedback. 



3.2. Initial conditions 

Given the elongated appearance of both the high and 
low velocity clouds in Serpens, we hypothesised the 
velocity structure of the Serpens cluster to be due to 
the collision of two elongated clouds (~ 1 pc long). 
Cloud collisions have long been thought to be an im- 
portant process in the ISM (since Por t 1954). They are 
potential triggers fo r star formation (e.g. |Scoville et al. 
1986[ Vallee|1995 ), and a re frequently found in galac 



tic scale simulations (e.g. |Dobbs|[2008[ |Tasker & Tan 
[2009 ). We do not suggest that cloud collisions are re- 
sponsible for all star formation in the Galaxy, but aim 
to test how well this scenario can model the specific 
case of Serpens. Taking individual clouds from galactic 
simulations to model the collision is beyond the scope 
of this paper, and furthermore does not allow any free- 
dom of the geometry of the clouds or collision. Instead, 
for our models we assume a simplified initial configura- 
tion with cylindrical clouds. The choice of cylinders is 
reasonable given that the most commonly observed as- 
pect ratio of molecular clouds is ~ 2 ( Koda et al. 2006 ). 
Thus our model is a simplified but plausible scenario in 
the local Galaxy. 

For our calculations, we based the properties of the 
cylinders on the observations of Serpens, but with the 
requirement that the cylinders are not too far from virial 
equilibrium. We required the star formation to be pri- 
marily driven by the collision, and preferably the cylin- 
ders should retain their elongated shape as much as pos- 
sible prior to the collision. Based on the observed spa- 
tial distribution of the Serpens' high and low velocity 
clouds, all the simulations we run involve two cylindri- 
cal clouds with radii of 0.25 pc. For all the calculations 
we have kept one cylinder vertical and the other tilted 
and, with the exception of run D (Table [T]), the cylinders 
have a length of 0.75 and 1 pc respectively. 

To ensure the stability of the cylinders prior to the 
collision, we used the formula for stability of finite 
cylinders given in |Bastien] ( |1983| ) to estimate the masses 
of the cylinders. 



GMf(L/D) 
LRgTlii 



(1) 



where Rg is the gas constant, G is the gravitational con- 
stant, M is the mass of the cylinder, L the length, T the 
temperature and jd the mean molecular weight (we as- 
sumed [I = 2.0 for these calculations). f(L/D) is a func- 
tion measuring the shape of the cloud, but is of order 
unity for the dimensions of our cylinders. The cylin- 
der is stable providing / < 0.8, a condition which 
reduces to approximately M(Mo)/T(K) < 0.8 for a 
1 pc cylinder. In addition, we include an external pres- 
sure term, which is subtracted from each SPH particle 
when calculating pressure gradients ( [Dale et al.|2009| ). 
This represents a confining pressure from a low density 
medium, which minimises diffxision of the outer parts 
of the cylinders, particularly in the runs including tur- 
bulence. 

To mimic the supra-thermal line widths observed 
in Serpens each model was also run for the same ini- 
tial configuration adopting a turbulent velocity field for 
the particles. The turbulent field was set up using the 
method of |Bate et al.l ( 2002 l, and described in depth in 
Dubinski et al . | ( T995| . The energy of the turbulent field 
E(k) is chosen to follow a power law of E(k) oc k~^, with 
k being the wavenumber. This gives, for a given scale, 
r, a velocity dispersion, cTy, such that CTy oc r^-^. This 
is similar to the observed Larson relations for clouds 
( Larson||1981 ), although we do not aim to capture the 
precise details of the internal motion, but rather to com- 
pare results with and without turbulence. 

From the velocity- space distribution of C^^ Q emis- 
sion in the Serpens Main Cluster (see Fig. [3 and'Duarte-l 
|Cabral et al.||2010l), we find that the higher velocity 
component (red in Fig. [Sj is tilted at 55° angle from 
the horizontal plane, and is mostly to the East on the 
map. The lower velocity component (blue in Fig. [3]) is 
roughly vertical and on the West of the map (overlap- 
ping only along the SE end of the filament of star for- 
mation). To match the observations, the higher velocity 
component corresponds to the longer and tilted cylin- 
der moving away from us on the right-hand side (see 
Fig. |4] red cylinder). This cylinder has an elevation an- 
gle of 55° above the horizontal plane and an azimuth 
angle of 0°, i.e. the axis is in a plane parallel to the 
Yz' plane. The lower velocity component corresponds 
to the vertical cylinder which is coming toward us from 
the left-hand side (see Fig.[4]blue cylinder). 

The velocities of the cylinders are chosen to agree 
with the velocity diff'erence of the gas in the Serpens 
cluster (see lower panel on Fig. [2]), which would cor- 
respond to the velocities of the clouds during the col- 
lision. The initial velocities of each cylinder are v = 
±1 kms"^ i.e. Vx' and Vy = ±0.77 kms"^ The cylinders 
have no initial net velocity in the z direction (and with- 
out turbulence the particles have zero velocity in the z 
direction). 

Each of the calculations were run with and with- 
out a turbulent velocity field. The diff'erences between 
the four runs presented in this paper are sketched in 
Figure |4j and summarised in Table [T] where "Collision" 
specifies if the collision is off'-centered or head-on with 
respect to the centre of the cylinders; "LOS < V" rep- 
resents the angle between the line of sight (LOS) and 
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Model A 



Top view 




Model B 



Top view 




Top view 



Model C 




Model D 



Top view 




Fig. 4. Illustrative diagrams showing the starting configuration of the cylinders for the direct collision (model A, top left panel), 
off'set collision (model B, top right panel; and model C, lower left panel) and off'set collision with longer cylinders (model D, 
lower right panel). The black dashed line represents the line of sight, and the green dashed lines are the direction of the motion 
of the cylinders. On the top right corner of each panel there is a sketch with the configuration as viewed from the top. 



the direction of the cylinders' motion; "Length" is the 
height of the longer (i.e. the tilted) cylinder in each 
run; "Velocities" specify the initial distribution of ve- 
locities: "Non-Turb", means zero velocity dispersion in 
the initial conditions and "Turb" indicates a turbulent- 
generated velocity dispersion amplitude of 0.5 kms"^ 
(in each direction). This is similar to the observed tur- 
bulent velocity dispersion in Serpens of the order of 
0.5 kms" HDuarte-Cabral et al.,2010j 

All other properties of the cylinders, such as the 
temperatures and densities, are fixed. In the "short- 
cylinder" runs (A, B and C) we used 500,000 par- 
ticles in total, with 250,000 particles in each cylin- 
der. Initially, the particles are uniformly randomly dis- 
tributed, to provide a roughly uniform density distribu- 
tion within each cylinder. 



As discussed in [Duarte-Cabral et aL] (120 101), the 
NW sub-cluster tends to have a uniform temperature of 
10 - 13 K. On the other hand the temperature of the SE 
sub-cluster is more difl&cult to estimate, given that there 
are tw o clumps along the line of sight. Duart e-Cabral| 
et al.j ( ^20 10^ obtained gas temperatures between 10 
and 20 K. However, as shown by Eq. [T] a kinetic 
temperature of 30 K is required for a 30 M© cylinder 
to be relatively stable. Hence we decided to add in the 
simulations this extra internal pressure in order to avoid 
the collapse of the cylinders before the collision. We 
chose masses of 30 and 45 M© for the two cylinders, 
and verified that for a temperature of 30 K the cylinders 
did not collapse prior to the collision. This does not 
have any impact on the velocity and column density 
structure discussed in the rest of the paper. 
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Table 1. Configuration of the diff'erent simulations. 



Run 


ColHsion 


LOS < V 


Length 


Velocities 


Anon- 


-T 


Direct 


45° 


Ipc 


Non-Turb 


At 




Direct 


45° 


Ipc 


Turb 


-^non- 


-T 


Off'set 


45° 


Ipc 


Non-Turb 






Off'set 


45° 


Ipc 


Turb 


Ct 


Off'set 


0° 


Ipc 


Turb 


Dnon- 


-T 


Off'set 


0° 


1.5pc 


Non-Turb 



We set the cylinders to collide at an angle of 45° rel- 
ative to the line of sight. In addition, we chose to posi- 
tion the cylinders so that either they collide head on (run 
A), or their major axis are off'set by 0.25 pc (run B). In 
the direct collision, model A, we centre one cylinder at 
cartesian coordinates (0,0,0) pc and the second (longer) 
cylinder at (0.8,0.8,0) pc (Fig.|4]top left panel). For the 
off'set collision, model B, the second cylinder is instead 
placed at (0.8,1.05,0) pc. We illustrate the off'set config- 
uration in Fig. |4] (top right panel). 

The results from the runs A and B are summarised 
in Sections [4. 1 . 1 1 to [4. 1 .2| These runs, however, show 
that by the time sink particles are formed the two cylin- 
ders fully overlap along the line of sight, even though 
they did not interact everywhere. This means that we 
can see two velocity components co-existing in most 
regions (e.g. Fig.[8]top left panel), inconsistent with the 
observations where the two velocity components over- 
lap only where the collision takes place. 

Therefore, we changed the perspective so that the 
motion of the cylinders is purely along the line of sight, 
model C (Fig. [4] bottom left panel), to restrict the re- 
gion where the two-components overlap to where the 
cylinders interact directly. Finally, model D (Fig.|4]bot- 
tom right panel) aims to reproduce the whole extent of 
the observed cloud and, therefore, it is similar to model 
C, except the cylinders were made one and half times 
longer by extending both the top and bottom of the 
cylinders. The number of particles was correspondingly 
increased in this calculation. 



4. Results 

As described in Section |2] the Serpens Main cluster 
is divided into two sub-clusters, the SE and NW sub- 
clusters. All the observational evidence tends to show 
that the SE sub-cluster is at the interface of a cloud- 
cloud collision, as opposed to the much more quiescent 
NW sub-cluster. For this reason, in the cloud-cloud col- 
lision simulations we present here, we first attempt to 
reproduce the physical properties of the SE sub-cluster. 
However, given the obvious physical connection be- 
tween these two sub-clusters we also try, in a second 
set of simulations, to trigger the formation of a NW-like 
sub-cluster, as a by-product of the collision. 

For comparison of the velocity structure of the sim- 
ulations w e will use the o bservational data and results 
from Dua rte-Cabral et al.| (2010): IRAM 30m telescope 
observations of C^^O J=1^0. We also refer to the anal- 



ysis of Serpens molecular data from the JCMT GBS 
HARP data of C^^O J=3^2 ( [Graves et al.||2010| ). For 
H2 column density comparisons, we used the SCUBA 
850 Jim emission from JCMT (Davis et al.|T999l ). 

The comparison of the simulations with the obser- 
vations is primarily focused on observational character- 
istics of the Serpens SE sub-cluster given that this is 
the region which appears to be directly influenced by 
the collision. However, even in the short-cylinder runs 
which do not form the NW sub-cluster, some of the ve- 
locity characteristics of the entire cloud are also taken 
into account. Thus, the following characteristics are the 
main focus points for comparison: 

- Column density structure: an elongated/filamentary 
shape aligned in a NW-SE direction, sub-clumped 
(Fig.[T]). 

- The mass and size of the southern clump compared 
with those of the denser parts of the SE sub-cluster 
(~30 M© within -0.025 pc^, measured using the 
850 fim emission above 190 mJybeam"\ which 
corresponds to column densities above 0.1 gcm"^ 
at a temperature of lOK). 

- Overall velocity structure: a single red- shifted com- 
ponent in the north, with a fainter blue component 
spatially off'set to east (top panel of Fig. |2] top pan- 
els and Fig. [3]). Double (overlapping) components 
on the south, where the sources are being formed 
(Fig. |2] lower panel and Fig. 3]) with a gradient from 
east to west of increasing velocities, up to 1.5 kms"^ 
(Fig. [13]). Moving further south, this velocity struc- 
ture should evolve back into one component along 
the less dense gas of the filament (Fig. 13 and lower 
PV diagram of Fig. 17 of [Graves et al.|p010| )). 

- Velocity structure of less dense material: thin east- 
ern filaments perpendicular to the main filament of 
Serpens (seen in blue and green on the left-hand 
side of the main filament in Fig. 13 ). 



In order to compare with the observations, we first 
constructed a datacube from the simulated 3D cloud 
collisions. For each simulation, we define the x direc- 
tion as the line of sight, so that the plane of the sky is 
the yz plane. Choosing the time- step of the simulations 
where the first sink particle is formed, we created a dat- 
acube of column density for a space- space- velocity 3D 
grid. In the spatial planes we convolved the models with 
a Gaussian of FWHM of 10 pix (0.02 pc) which corre- 
sponds to the 22" spatial resolution of the IRAM-30m 
telescope observations of C^^O J=1^0 at the distance 
of Serpens. While to reproduce the thermal velocity dis- 
persion of Serpens, the velocity space was convolved 
with a normalised Gaussian of 0.4 kms"^ full width half 
maximum (FWHM), correspondent to the thermal line 
width of H2 at lOK. 

Gas temperature, density and abundance all aff'ect 
the relationship between the true column density of 
a cloud and the emission seen in a molecular line. 
However as discussed in ,Duarte- Cabral et al.l ((2010]) 
and Grave s et al.| ( |201Ql ) in Serpens the C^^O appears 
to be a faithful tracer of the overall velocity structure of 
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Fig. 5. Colour scale and contours of the total column density along the line of sight, for the non-turbulent runs: the centered 
collision, Anon-x, on the left; and off-centered collision, Bnon-T, on the right. The contours levels are 0.01, 0.03, 0.1, 0.3 and 0.5 
g cm~^ for both figures. The white horizontal lines show the cuts for the PV diagrams in Fig.jS] 
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Fig. 6. Position- velocity diagram at z = -0.02pc (colour scale and contours) for the non-turbulent direct and off'set runs, Anon-x 
and Bnon-T, left and right respectively. The displayed velocity corresponds to the component along the line of sight. Contours at 
10~^, 10""^, 5 X 10~^, 10~^, 5xl0~^, 10~^ and 3 xlO~^ gcm~^. The red crosses represent the column density weighted velocities, 
and are plotted as an auxiliary tool to see the velocity changes along the diagram. Note that the velocities from the direct 
collision (left) in the central region seem to be more complex than those observed (Fig.[2]lower panel), whereas those from the 
off'set collision (right) show a smooth shift in velocities from one component to another, likely because in the off'set collision 
we see more clearly the parts of the clouds which do not collide. 



the cloud and not significantly afifected by outflow or in- 
fall motions, while globally the 850 jim emission traces 
the mass distribution in the region. To assess the success 
of the simulations in modeling Serpens, we therefore 
compare the models with the velocity structure of the 
C^^O and the overall column density distribution from 
the dust emission. 



4.1. Modeling the Serpens SE sub-cluster 

We initially try to reproduce the characteristics of the 
SE sub-cluster, using the short-cylinders calculations: 
runs A, B and C. In this section, we describe the re- 
sults of these simulations. Taking the best-fit configu- 
ration for the short-cylinders, we then perform simula- 
tions with more extended cylinders designed to repro- 
duce the whole cloud (Section 142]). 



4.1 .1 . Non-Turbulent Runs: Anon-x and Bnon-x 

Fundamentally, the only diflference between these two 
simulations is that in the Anon-T case the centres of 
gravity of the cylinders are colliding head-on, maximis- 



ing the loss of kinetic energy in the shock, while in the 
other case, Bnon-T, the centres of gravity are shifted and 
therefore the collision becomes softer, with the gas be- 
ing able to conserve more of its initial kinetic energy. 
While the first configuration is a very peculiar one, the 
second one is very general, corresponding to the most 
likely situations of cloud-cloud collision. As a result, 
the formation of sink particles in Bnon-T is slower than 
in Anon-T, and while it occurs only at the interface re- 
gion for the latter, in Bnon-T sink particles are more 
widely distributed, along a dense filam ent c onnecting 
both clouds (see Fig.|5]and App.[A] Fig. A.l ). 



The most significant diflTerence between the results 
of these two simulations is the velocity structure at 
the collision interface. In Anon-T, three distinct veloc- 
ity components are present: the two velocities of the 
individual cylinder, plus an intermediate velocity of the 
shocked interface gas (Fig. [6] left panel). In Bnon-T, the 
velocity structure is continuously connecting the initial 
velocities of the cylinders (Fig.[6]right panel). Thus, the 
offset collision (B) is better able to reproduce the obser- 
vations (Fig. [2]). 
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Fig. 7. Colour scale and contours representing the total column density along the line of sight, for the turbulent runs: the 
centered collision, At, on the left; and off-centered collision, Bt, on the right. The contour levels are as in Fig.|5] The white 
horizontal lines show the cuts for the PV diagrams shown in Fig. [8] Both turbulent runs show a final cloud which is much more 
structured better fitting the observations than Fig. [5] 
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Fig. 8. Position- velocity diagrams in colour scale and contours, at z = 0.16 pc (top), z = -0.06 pc (middle) and z = -0.30 pc 
(bottom) for the turbulent direct and off'set runs. At and Bt, left and right respectively. The contour levels and red crosses are as 
in Fig.|6] As expected from Fig.|6] the velocities are more complex in the A run than those observed. The top panels show that 
the direct collision (left) produces two velocity components in the northern part of the cloud, whereas both the off'set collision 
B (right) and the observations only show one (Fig. [2]). From the middle and lower panels, the diff'erences are less obvious, but 
the middle panel shows that the direct collision shows a more "messy" velocity structure rather than a clear transition from 
high to low velocities. 

4.1.2. Turbulent Runs: At and Bt 

In these simulations the net effect of turbulence is to 
render the cylinders inhomogeneous both in density and 
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velocity, making them more realistic looking clouds. 
One of the important feature we see in Fig. [7] (and 
App.[A| , Fig. |A.2| ), is the presence of channels of mate- 
rial, perpendicular to the main filament. These are rem- 
inisc ent of feature s often observed in molecular clouds 
(e.g. |Myers|[2Q09] ). These features will be further dis- 
cussed in Sec.|5J Another important aspect of injecting 
turbulence is that is generates density seeds. Before the 
first sink particles start to form, we see several regions 
where the density is starting to grow. As soon as one of 
these regions starts to be dense enough, it will attract 
the remaining less dense sub-clumps to form a major 
single structure. 

At the end of both runs, the material at a column 
density greater than 0. 1 gcm"^ is distributed over simi- 
lar sized regions to the non-turbulent models, 0.025 pc^ 
and 0.032 pc^ for runs At and Bj respectively. For 
At, a total of 29 M© is in these high column density 
regions (39% of the total mass of the cylinders), while 
for Bt the mass is 34 M© (45% of the total mass). 

In terms of the velocity structure, the direct colli- 
sion model. At, still produces too many distinct veloc- 
ity components with no smooth trends (Fig. [8] left pan- 
els), for the same reasons as for the Anon-x model. On 
the other hand, the Bt model produces a very interest- 
ing velocity structure. To the north (positive z) we see 
mainly the tilted cylinder and we begin to detect the 
vertical cylinder on the left-hand side. The transition 
between the two is not sharp nor double peaked, but a 
rather broad smooth transition (Fig. [8] top right panel). 
A well defined double peak structure is only detected 
in the central region of the model where sink particles 
have formed (Fig. [8] middle right panel). South of this 
region, the column density becomes again dominated 
by one component from the tilted cylinder traveling 
away from us (Fig. [8] bottom right panel), as also found 
towards the s outh end of the Serpens filament by'Graves 
|et al.| ([20T0|. The results from this run (Bt) are more 
consistent with the observations of the SE sub-cluster 
of Serpens (Sec. [2]), but the change in velocity across 
the cloud is smaller than that observed (-0.5 kms"^ in- 
stead of ~1 kms"^). 

We have also run a model with an intermediate 
turbulence level (not shown), with a velocity disper- 
sion of 0.3 kms"^. This model showed that the geo- 
metric configuration is the dominant factor in deter- 
mining the main velocity characteristics of the resulting 
cloud, while the density distribution is sensitive to the 
level of turbulence. Lower levels of turbulence result 
into a more filamentary structure, with more sink par- 
ticles forming along the denser filament. Increased val- 
ues of initial turbulence tend to disrupt the distribution 
of cloud's material, inhibiting the formation of further 
clumps. 



4.1 .3. Line of sight, short cylinder Run: Ct 

Model Ct is simply a re-projection of model Bt so that 
it is viewed along the axis of the relative motion of the 
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Fig. 9. The total column density for the purely along the line 
of sight turbulent run, Ct- Unlike runs A and B, we find a low 
density region in the north (large z), as with this orientation 
we see regions which have not collided. The contour levels 
are as in Fig.|5] 



cylinders. The time evolution of Ct is shown in App.[A| 



Fig. A.3 



From this perspective (Fig. [9]) the low density mate- 
rial is less filamentary and more spatially extended, as 
the parts of the cylinders which do not interact do not 
overlap along the line of sight any longer. Therefore, 
the extent of the lower density material is greater to 
the north, where the cylinders do not merge, than in 
the south (negative z). In this projection, the material at 
high column density above 0. 1 gcm"^ covers a region 
of 0.02pc^ in area with a mass of 22M0. 

The general shape and trend of the velocity struc- 
ture (Fig. [To]), although slightly more complex, is sim- 
ilar to that seen in Bt and consistent with the velocity 
trend seen on the observational PV diagrams. Overall, 
this run (C) reproduces Serpens better than run B, as it 
now shows a more evident velocity change across the 
cloud. 



4.2. Representation of both SE and NW Serpens 
suh-clusters 

4.2.1 . Line of sight, long cylinders, Run: Dnon-T 

Overall, model Ct satisfactorily reproduces the obser- 
vations of the SE sub-cluster. However it does not have 
enough Jeans masses to form a second sub-cluster com- 
parable to the NW sub-cluster in Serpens even though 
the velocity structure of the simulation is very similar to 
that observed. Therefore, we increased the length and 
mass of the cylinders to investigate the possible for- 
mation of a separate structure in the north where the 
cylinders do not interact. If the non-interacting north- 
ern region is massive enough to be close to unstable, 
the collision in the south of the cylinders may induce 
its rapid collapse without greatly aff'ecting its systemic 
velocity. 

A non-turbulent model with -1.5 pc long cylinders, 
colliding off'-centre and along the light of sight (model 
Dnon-T ) does indeed form two sub-clusters. The colli- 
sion in the south leads to the formation of a sink particle 
in the south, i.e. the synthetic SE sub-cluster. Later on. 
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a sink particle is formed in the north of the tilted cylin- 
der, along the cylinder axis (time evolution in App. [A| 
Fig. |A.4| ). Finally, more sink particles are formed in the 
south, as the vertical cylinder crosses through the tilted 
one (Fig. [TT]). At the end of this run, when sink parti- 
cles are formed throughout the clouds, the column den- 
sity distribution is much more filamentary than for the 
short cylinder calculation. There are two visible sub- 
clusters: one in the NW and one in the SE. The total 
mass in these sub-clusters amounts to 57 M©, 52% of 
the total mass in this case, distributed over an area of 
0.035 pc^. However, the relative size and mass of each 
individual sub-cluster is not quite as similar as seen in 
Serpens. The simulated SE sub-cluster has a mass of 
about 30 M© in an area of -0.015 pc^, whereas the NW 
sub-cluster contains 12 M© in a similar area. 

The formation of a sink particle in the north occurs 
at 8 X 10^ yr and is mainly due to the initial instability 
of the cylinder itself. As a test, we ran a simulation with 
only an isolated non-turbulent cylinder, with the same 
mass and size as the tilted one. This simulation showed 
that it would form sink particles along the cylinder axis 
at lO^yr. Therefore, the collision in the south does not 
trigger the formation of this sink particle in the north, 
but only speeds up its collapse. 

Interestingly, several diff'erent turbulent runs with 
this configuration failed to induce the collapse of any 
structure in the north. The collapse only occurs with 
extremely low levels of turbulence (with amplitude of 
~ 0.05 kms"^), almost indistinguishable from the non- 
turbulent run. This is because the turbulent velocities 
cause the cloud to disperse from its initial configura- 
tion and become less gravitationally bound than in the 
absence of turbulence. From this, we conclude that the 
eff'ect of the collision felt by the north region is quite 
subtle, and would only have an influence on the col- 
lapse there if the region is almost gravitationally unsta- 
ble prior to the efl'ect of any external perturbation. 

The discrepancy in the relative masses of the two 
sub-clusters in the simulation compared to Serpens, and 
indeed the failure to form two sub-clusters in the tur- 
bulent simulations, could both be addressed by relax- 
ing the over-simplistic uniform conditions in our initial 
conditions. Both the density distribution and turbulent 
velocity distributions are likely to be more inhomoge- 
neous in real colliding clouds than in our models. Also, 
this discrepancy could be suggesting that we are miss- 
ing a physical ingredient such as magnetic pressure, 
which could help maintaining the filamentary shape of 
the clouds even with high level of turbulence. However 
we refrained from adding such additional complexity 
to the models as the nature of the inhomogeneities is 
poorly constrained and unlikely to provide any further 
significant physical insight into the processes in this re- 
gion. 

In terms of its spatial evolution, the velocity struc- 
ture of the Dnon-T is similar to the observed one. From 
Fig. 12 we see that as we move from north to south 
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(i.e. from positive z to negative z), two spatially sepa- 
rated components, which then overlap where the col- 
lision takes place, to finally end-up with a single ve- 



Fig. 10. Position- velocity diagram in colour scale and con- 
tours, at z = 0.14 pc (top), z = -0.07 pc (middle) and 
z = -0.30 pc (bottom), for the turbulent, line of sight, short 
cylinders run Ct- The contour levels and red crosses are as in 
Fig. [6] Compared to model Bt (Fig.[8]right column), we see 
more evidence of the original velocities of the clouds from the 
non-interacting parts, due to the change of perspective. This 
calculation produces the best fit to the observations (Fig. [2]). 



locity component in the south. These PV diagrams can 
be compared to the m ore detailed PV diagrams ac ross 
the entire cloud from [Duarte-Cabral et aL] ( |2010| ) and 
[Graves et al.| ( |2010j) . Also note that the column densi- 
ties are higher at velocities corresponding to the tilted 
cylinder. An increase in the density of the tilted cylin- 
der to produce a more massive NW sub-cluster would 
increase the diff'erence in column densities between the 
two velocity components, more closely matching the 
observations. 



4.3. Timescales 

One interesting question for the star formation history 
of Serpens is to confirm whether or not the older popu- 
lation of pre-main sequence star of 2 Myr of age could 
be formed in the same cloud-cloud collision event. For 
this we need to estimate a number of characteristic 
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Fig. 11. Total column density for the non-turbulent, purely 
along the line of sight run with longer cylinders, Dnon-x, to 
test the production of both the NW and SE sub-clusters of 
Serpens. At this stage, sink particles have indeed formed 
where the collision is happening and in the north-west (within 
the tilted/longer cylinder). Given the geometry of the colli- 
sion, the two cylinders do not end up colliding in the north. 
However, the only perspective where we do not see the cylin- 
ders overlap is the one we took for this run, where the line of 
sight is aligned with the motion of the cylinders. The contour 
levels are as in Fig.|5] 



timescales. First, all the simulations correspond to a to- 
tal elapsed time of ~ 8 - 9 x 10^ years. However, it 
is only after about 2 - 3 x 10^ years, about a third of 
the total simulation time, that the two cylinders start to 
interact. This time delay is needed, especially in the tur- 
bulent cases, to allow the cylinders to relax their initial 
density distributions. The time from the start of the in- 
teraction between the two clouds until they start form- 
ing sink particles is of the order of 6 x 10^ years. For 
the 2 Myr stellar population to be able to form through 
this collision, the duration of the interaction time has to 
be at least 2 Myr. In fact, we would need the cloud to 
be an order of magnitude larger (i.e a few parsec wide) 
in order to still have an ongoing collision 2 Myr after 
the formation of the first protostars. This seems unre- 
alistic, rather the cloud collision scenario we present 
for Serpens is consistent with the idea of two separate 
bursts of star formation in the region. 



5. Discussion 

Overall the simulation which best represents the SE 
sub-cluster of Serpens star forming region is the model 
Ct, the off'set turbulent model with short cylinders, re- 
projected so that the line of sight is coincident with 
the direction of motion. It has a velocity structure sim- 
ilar to that observed, both comparing PV diagrams and 
the general trend on the velocity coded 3 -colour plots 
(Fig. [13] and Fig. [14]). Note that this simulation only 
represents the southern part of Serpens as the NW 
sub-cluster is not directly involved in the collision (cf 
Sec. 




4.1). Despite this, the 3-colour plot of this run 



(Fig.|l4[) shows a good agreement with that of the obser- 
vations (Fig. 13). Model Ct shows an overall velocity 
gradient of ~2 kms"^ over 0.2 pc"^ from blue in east- 
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Fig. 12. PV diagrams of the long cylinder run Dnon-x at z = 
0.16 pc, z - -0.02 pc, z - -0.09 pc, z = -0.21 pc and z - 
-0.37 pc (from top to bottom). Contour levels and red crosses 
are as in Fig. [6] Despite the lack of structure typical from 
a non-turbulent case, these agree with the more detailed PV 
dia grams from observat ions from |Duarte-Cabral et al.| ( [20T0| 
and |Graves etari ( [20T0l ). 
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Fig. 13. Velocity coded 3 colour plot of Serpens from 
the C^^O J=3^2 data from the JCMT GBS. Each colour 
represents the maximum value in the velocity intervals: 
blue: 5 1.1 kms"^; green: 1-1 ^ 8.3 kms~|; and red: 
8.3 ^ 11 kms~^ Also discussed in Graves et al. (2010i. The 



two sub-clusters of Serpens are indicated by the black el- 
lipses. 
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Fig. 14. Velocity coded 3 colour plot of the simulation Ct- 
Each colour represents the total column density summed over 
the following line of sight velocities: blue: 0.3 2 kms~^; 
green: -0.3 0.3 kms~^; and red: -2 -0.3 kms~^ This fig- 
ure is to be compared with Fig. [13] particularly around the 
SE sub-cluster. The synthetic SE sub-cluster is indicated by a 
black ellipse. 



ern regions to red on the west (Fig. 14 ), similar to what 
is observed. 

The 3 -colour velocity figures also show the resem- 
blance of the simulation Ct with the observations in the 
way that the filament extending south is mostly repre- 
sented by the red/green velocities, while the blue part 
is mainly seen on the denser parts where the stars are 
being formed. The green-blue component is also seen 
further east, forming some less dense filaments perpen- 
dicular to the main filament. Note that we do not see 
these types of filament on the red side. 

Understanding the nature of these filaments can be 
important. If these filaments had been caused by the 
initial turbulence, then we would expect these to exist 
on both directions. Alternatively, these filaments could 
be the consequence of some particles left behind as the 
cylinders move, resembling a tail. However, in simula- 
tion C, the cylinders' motion is along the line of sight, 
so if these filaments resulted from a tail of material left 
behind during the collision they would be behind or in 
front of the main structure. Their appearance therefore 
indicates that these perpendicular filaments are caused 
by the geometry of the collision, in particular, its asym- 
metric nature. These filaments place a strong constraint 
to the geometry of the collision. The geometry of simu- 
lation A (as well as other tests not shown) do not repro- 
duce them. 

In terms of its total column density distribution, 
22 M© in 0.02 pc^, model C closely matches the SE 
sub-cluster. Even the projected distribution of material 
above a column density of 0.1 gcm"^ resembles the 
shape of the SE sub-cluster. 

To start to reproduce the NW sub-cluster in addition 
to the SE sub-cluster requires the longer non-turbulent 



clouds used in model Dnon-T- In this model a region 
close to being gravitationally unstable can be perturbed 
and its collapse hastened by the collision, even though 
it is not directly involved in the collision. This second 
sub-cluster forms with the velocity of its native cloud 
and collapses smoothly and independently of the 
southern region. If this model is allowed to run further, 
the NW sub-cluster falls onto the SE sub-cluster. 
Intriguingly, recent observations of the magnetic field 
in Serpens ( [Sugitani et al.|[2QTQ] ) appear to suggest 
the start of such a collapse which could result in the 
merging of the two sub-clusters. 



6. Conclusions 

Serpens is a very interesting star forming region, not 
only for its youth, but also for the striking diflTerences 
between the two sub-clusters that compose the active 
star-forming portion of the cloud. Even though they 
are at similar stages of evolution, with most sources 
between Class and Class I protostars and similar 
dust continuum properties, the gas emission reveals that 
these have not only diflTerent velocity characteri stics but 
also diflTerent temperature distributions (Duarte-Cabral 
let al.|2010| ). 

Motivated by the two velocity components seen in 
the southern sub-cluster of Serpens, and the higher tem- 
peratures detected there, we performed several SPH cal- 
culations of cloud-cloud collisions. The configurations 
used for these simulations were based on the observed 
morphologies. 

For the SE sub-cluster, a model of two colliding 
clouds is able to reproduce both the column density 
structure, a centrally condensated filament aligned in 
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a NW-SE direction, and the two velocity components 
seen where the star formation driven by the colHsion 
is occurring. The same simulations did not, however, 
produce a second sub-cluster, similar to the NW sub- 
cluster of Serpens. Therefore, this sub-cluster does not 
seem to be the direct result of the collision. This was 
already suggested by the NW sub-cluster's well "be- 
haved" temperature profile and velocity structure, as 
well as the uniform age of sources within the sub- 
cluster. However, the similar stage of evolution of the 
sources from the two sub-clusters and their proximity, 
suggests that the two events are not totally independent. 

A simulation with more elongated cylinders and in- 
creased masses provides a possible explanation. The 
presence of a marginally stable region in the northern 
part of one of the colliding clouds can have its collapse 
induced and quickened by perturbations driven by the 
cloud-cloud collision. 

We consider a cloud-cloud collision scenario to be 
the best description of the driving of the star forma- 
tion history in Serpens. Not only can it reproduce the 
observed velocities and column densities, as it off'ers a 
plausible explanation for why the two sub-clusters are 
so similar in some aspects and yet so diff'erent in others. 
Although cloud rotation may produce similar general 
velocity gradients to those observed, the complexity of 
the region is better explained with such a collision sce- 
nario, which is in essence similar to a shear-motion also 



suggested by |Qlmi & Testi| ( |2002| ). 

Despite the successful scenario provided by a 
cloud-cloud collision model, we failed to reproduce 
all of the Serpens characteristics in one single run. 
Additional support against gravity is required in order 
to sustain the existence of two diff'erent sub-clusters as 
in Serpens. The existing magnetic field of the region 
( ,Sugitani et al.,2QlQj could provide such a support. 
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Appendix A: Time snapshots of tlie models presented 



Figures [AJ] to | A.4| show three snapshots to illustrate the time evolution of the models presented in this paper. Each 
frame is a projected column density plot as seen from the chosen line of sight. The time (in years) is show in the 
top-right corner of each frame. Figure A.l shows three frames for the two non-turbulent runs presented, Anon-T 
and Bnon-T- Figure A.2 shows the time evolution for the same models, A and B, but with turbulence. Figure A.3 
shows the snapshots of model C, which is essentially the same as model Bj only with a different line of sight. 
For these three runs, the last frame corresponds to the formation of the first sink particle and is the frame used 
to construct a datacube for comparison with the observations. Finally, Fig. A.4 shows three snapshots for run D. 
Note that for this run we are showing frames after the formation of the first sink particle, to demonstrate where 
and when the second clump in the north is able to form a sink particle. In this case, this last frame was the one 
used for the comparison with the observations. 




I 



Fig. A.l. Three time snapshots of the total column density along the line of sight, for the non-turbulent runs: the centered col- 
lision, Anon-T (top) and off-centered collision, Bnon-T (bottom). The three frames correspond to the beginning of the simulation 
(first frame), as soon as the cylinders start to collide (second frame) and when the first sink particle is formed (last frame). 
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Fig. A.2. Three time snapshots of the total column density along the line of sight, for the turbulent runs: the centered collision, 
At (top) and off-centered collision, Bt (bottom). The three frames correspond to the beginning of the simulation (first frame), 
as soon as the cylinders start to collide (second frame) and when the first sink particle is formed (last frame). 
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Fig. A.3. Three time snapshots of the total column density along the line of sight for Ct- The three frames are as before. Note 
that this is the exact same model as Bt, but viewed from a diff'erent perspective. 
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Fig. A.4. Three time snapshots of the total column density along the line of sight for Dnon-T- The three frames are at the start 
of the simulation (left), when the first sink particle forms in the south (middle) and when a sink particle is formed in the north 
(right). 
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